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I. INTRODUCTION. 



BES Collaboration observed a significant threshold enhancement of pp mass spectrum in the radiative decay J 
jpp [ij. Recently BES CoUaboration reported the results on X{1835) in the J/ip ^ -f rj'n+TT' channel Under the 

strong assumption that the pp threshold enhancement and X{1835) are the same resonance, Zhu and Gao suggested 
X(1835) could be a J^^ = 0"+ = 0+ pp baryonium 

Theoretical work speculated many possibilities for the enhancement such as the t-channel pion exchange, some kind 
of threshold kinematical effects, as new resonance below threshold or pp bound state [sl-fl^. 

In a series of papers (isl - llTj a method has been developed which is convenient for analysing relativistic three-hadron 
systems. The physics of the three-hadron system can be described by means of a pair interaction between the particles. 
There are three isobar channels, each of which consists of a two-particle isobar and the third particle. The presence 
of the isobar representation together with the condition of unitarity in the pair energies and of analyticity leads to a 
system of integral equations in a single variable. Their solution makes it possible to describe the interaction of the 
produced particles in three-hadron systems. 

In our papers relativistic generalization of the three-body Faddeev equations was obtained in the form of 

dispersion relations in the pair energy of two interacting quarks. The mass spectrum of S'-wave baryons including u, d, 
s quarks was calculated by a method based on isolating the leading singularities in the amplitude. We searched for the 
approximate solution of integral three-quark equations by taking into account two-particle and triangle singularities, 
all the weaker ones being neglected. If we considered such approximation, which corresponds to taking into account 
two-body and triangle singularities, and defined all the smooth functions of the subenergy variables (as compared 
with the singular part of the amplitude) in the middle point of the physical region of Dalitz-plot, then the problem 
was reduced to the one of solving a system of simple algebraic equations. 

In the previous paper [2l[ the relativistic six-quark equations are found in the framework of coupled-channel 
formalism. The dynamical mixing between the subamplitudes of hexaquark are considered. The six-quark amplitudes 
of dibaryons are calculated. The poles of these amplitudes determine the masses of dibaryons. We calculated the 
contribution of six-quark subamplitudes to the hexaquark amplitudes. 

In the present paper the relativistic six-quark equations including u, d quarks and antiquarks are found. The 
nonstrange barionia BB are constructed without the mixing of the quarks and antiquarks. The relativistic six-quark 
amplitudes of the baryonia are calculated. The poles of these amplitudes determine the masses of baryonia. In 
Sec. II the six-quark amplitudes of baryonia are constructed. The dynamical mixing between the subamplitudes of 
baryonia is considered. The relativistic six-quark equations are obtained in the form of the dispersion relations over 
the two-body subenergy. The approximate solutions of these equations using the method based on the extraction of 
leading singularities of the amplitude are obtained. Sec. Ill is devoted to the calculation results for the baryonia mass 
spectrum and the contributions of subamplitudes to the baryonia amplitude (Tables I, II, III, IV). In conclusion, the 
status of the considered model is discussed. 



•Electronic address: gerasyuta@SG6488 . s pb.edu 
t Electronic address: matskev@pobox.spbu.ru| 



2 



II. SIX-QUARK AMPLITUDES OF THE BARYONIA. 



As explained in the previous paper j2l| the relativistic generahzation of the three-body Faddeev equations was 
obtained in the form of dispersion relations in the pair energy of two interacting quarks. The pair quarks amplitudes 
qq qq are calculated in the framework of the dispersion N/D method with the input four-fermion interaction [23j 
with quantum numbers of the gluon [U, [2^ . 

The construction of the approximate solution is based on extraction of the leading singularities are close to the 
region Sik ~ 4m^. Such a classification of singularities makes it possible to search for an approximate solution of 
equations, taking into account a definite number of leading singularities and neglecting the weaker ones. 

We derive the relativistic six-quark equations in the framework of the dispersion relation technique. We use 
only planar diagrams; the other diagrams due to the rules of 1/Nc expansion [26l - l28l l are neglected. The current 
generates a six-quark system. The correct equations for the amplitude are obtained by taking into account all 
possible subamplitudes. It corresponds to the division of complete system into subsystems with a smaller number of 
particles. Then one should represent a six-particle amplitude as a sum of 15 subamplitudes: 



(1) 

This defines the division of the diagrams into groups according to the certain pair interaction of particles. The 
total amplitude can be represented graphically as a sum of diagrams. We need to consider only one group of diagrams 
and the amplitude corresponding to them, for example Ai2. We shall consider the derivation of the relativistic 
generalization of the Faddeev- Yakubovsky approach. 

In our case the low-lying baryonia are considered. We take into account the pairwise interaction of all quarks and 
antiquarks in the baryonia. 

The system of graphical equations Fig. 1 is determined using the selfconsistent method. The coefficients are 
determined by the permutation of quarks [H, Is^l ■ We should discuss the coefficient multiplying of the diagrams in 
the equations of Fig. 1. For example, we consider the first subamplitude A}"" (s, S12345, S1234, S123, S12). In the Fig. 1 
the first coefficient is equal to 2 (permutation particles 1 and 2). The second coefficient is equal to 6 = 2 (permutation 
particles 1 and 2) x3 (we consider the third, the fifth, the sixth particles). The similar approach allows us to take 
into account the coefficients in all equations. 

In order to represent the subamplitudes Al , A\ , A\ , ^ , ^ ^ in the form of dispersion relations, 
it is necessary to define the amplitudes of qq and qq interactions. 

We use the results of our relativistic quark model [25j and write down the pair quark amplitudes in the form: 



an{Sik) 



1 - Bnis^k) 



(2) 



B^s.,)^ f ^PAffm^A, (3) 

J TT S^^ - S,k 

In the case of the s-channel amplitudes we use the matrix element: (^q'^O^q'^ ^ (^q'^ O^q'^^ , where qc = qC is the 

charge-conjugated spinor. are operators of different types of the four-fermion interaction (i = S,V,T, A, P), a, 6, 
a', b' flavor indices. 

Here Gn{sik) are the diquark vertex functions (Table V). The vertex functions are determined by the contribution of 
the crossing channels. The vertex functions satisfy the Fierz relations. Since the vertex functions depend only slightly 
on energy it is possible to treat them as constants in our approximation and determine = N oi N/D method. 
These vertex functions are generated from gluon coupling constant g. Bn{sik) and Pn{sik) are the Chew-Mandelstam 
functions with cutoff A [31] and the phase spaces, respectively: 



Pn{s^k,J^^) = [a{n,.r^)- i^+/3(n,J^^)+<5(n,J^^) 

yj {sik - (mi + mkY){sik - {rrii - rukY) 
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The coefficients a(n, J^*^), f3{n, J^'~') and J^'~^) are given in Table V. 

Here n = 1 coresponds to gg-pairs with = 0^, n = 2 corresponds to the gg-pairs with = 1^, ri = 3 defines 
the gg-pairs with — 0+ , n = 4 coresponds to = 1+ gg- pairs. 

Let us extract two- and three-particle singularities in the amplitudes A\ (s, S12345, S1234, Si23, S12), 

^1 (s,_Si2345, ■S1234, S123, S12), ^1 (s, S12345j S1234, S123, S12), ^2 ^ (s, Sl2345j Sl234, Sl2, S34), 
A^""l" 1 (s,Si2345,Sl2, 534,355): 



.1"". X (s, S12345, S1234, S123, Sl2)Sl""(sl2) 

^1 S12345, S1234, S123, •Si2j = j- J5 -, , (O) 

, N (s, S12345, S1234, Sl235 Si2)i3id<i(si2) , r.\ 

^1 (s,Sl2345,Sl234,Sl23,Sl2j = r; -. , (6) 

.l"<f/ N Ol\ (s, S12345, ■S1234, S123, Sl2)-Bi«d(Si2) , , 

^1 (s,Sl2345,Sl234,Sl23,Sl2j = r; -, , (7) 

■ 1""1<I<J, ^ a2""^'^''('S, Sl2345,Sl234,Si2,S34)-Bi>.,.(si2)Bi<id(s34) 

A2 (S, Si2345,Sl234,Sl2,S34j = j-j 5 7 TTT^ g 7 7^ , (oj 



d. -1 dd 



^3 ^ ^ (s,Si2345,Sl2,S34,S56) 



03 (S, S12345, gl2, g34, gse)-^!^" (5l2)-Bi„d (534)-Bid-d-(s56) 

[1 - Bi..(si2)][l - Si„d(s34)][l - 5id-d(s56)] 

(9) 



We do not extract four-particles singularities, because they are weaker than two- and three-particle singularities. 

We used the classification of singularities, which was proposed in paper [32| . The construction of the approximate 
solution of Eqs. (5) ~ (9) is based on the extraction of the leading singularities of the amplitudes. The main 
singularities in sik — {nii + rrikf' are from pair rescattering of the particles i and k. First of all there are threshold 
square-root singularities. Also possible are pole singularities which correspond to the bound states. The diagrams of 
Fig. 1 apart from two-particle singularities have triangular singularities and the singularities defining the interactions 
of four, five and six particles. Such classification allows us to search the corresponding solution of equations by 
taking into account some definite number of leading singularities and neglecting all the weaker ones. We consider the 
approximation which defines two-particle, triangle and four-, five- and six-particle singularities. The contribution of 
two-particle and triangle singularities are more important, but we must take into account also the other singularities. 

The five functions ai are the smooth functions of Sik, Sijk, Sijki Sijkim as compared with the singular part of the 
amplitudes, hence they can be expanded in a series in the singularity point and only the first term of this series should 
be employed further. Using this classification, one defines the reduced amplitudes ai as well as the B-functions in the 
middle point of physical region of Dalitz-plot at the point sq: 
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^0 - , (10) 

i,k = l 



Sl23 



3 3 
So > mfi. - > TO? 



i<k 



i=l 



(11) 
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S1234 = So X! '"^fe - 2 XI (12) 



,fc=i 
i<k 



Such choice of point sq allows us to replace integral equations (Fig. 1) by the algebraic equations (13) - (17), 
respectively: 
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af" = A + 27i(l"«l"")af" +67i(l""l"'^)ar , (13) 
af = \ + 2h{l'^h'''^)af + &h{l'^h^'^)af , (14) 
af ' = A + 27i(l"<^"l"")af'' + 2h{Vh'^'^)af^ + 47i(l"'^"l«'^")af ' + 472(l"'^"l""l'^"'^")ar"''' , 

(15) 

api'" = A + 274(l""A"")ar" +274(l'^"'^>"l'^"'^")ar +473(1""1'^"'^>'^")^ 



+ 476(l""l'^"'^"l""l'*"'^")ar"^'' +478(l""l'^'^l""l"'^l'*'^)ar"^"''^''' . 
^^uu^u4^44 ^ 279(l«"l"'^"l^"^"l"")af " + 279(l<^'^"l«'^"l«''l'^"'^")a}'" + 279(l"''l«'^"l^"^"l"'^>f' 
+ 479(l""l'^''^"l"'^"l"'^")af' +279(l'^"'^"l"'^"l""l"'^")af' + 47io(l""l"^"l'^"'^"l""l'^"'^")ar"^' 
where Aj are the current constants. We used the functions 7i, I2, I3, I4, Iq, h, hi ho- 



2 

(mi+m2)^ —1 —1 



23(2) + 

X / dz3{2) ^ 



hiijk) 



- zf{2) - zl{2) - zl{2) + 2^1(2)^2(2)^3(2) 

Z3(2)- 

1 1 

l-i?,(sl3)l-i?fc(s^4)' 

(mi+n.2)^Ai 



(mi+mg)^ 

(m3 + m4)^Aj 

rfs^4Gi(«o')/'j(-<34) 7^^1(3) ydZ2(3) 



Is{ijklm) 



(mi+Tn2) Aj 



(mi+m2)^ 



(16) 



(17) 



''^'^> B,{sl^) J n s[,-sl' J 2 1-B,(5l3)' ^ ^ 

(mi+m2)^ —1 

ir^^ + r^ +1 +1 

_ B,(5;3)B,(g24) r ds[, GUsl')Pi{s[2) 1 Z" rf£l(2) I" d^2(2) 



(19) 



f ds'3^G^^{s-^-)p,{s'^^) f dzi{3) f dz2{3) 1 

J s'^^-sf J 2 j 2 1-Sfe(43)' ^ ^ 

(m3+m4)^ —1 —1 

h{ijk) = h{ik), (21) 
hiijkl) = h{ik)-h{jl), (22) 
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X 



(m3+r,.4)^Aj 



^34 '=0 



+ 1 +1 +1 Z4(8)+ +1 Z6(8) + 

1 [ dzi{S) fdz2{8) f dz3{8) f ^ f dz5{8) 



dzi{8) / / dz6{8) 



(27r)2 7 2 J 2 J 2 J J 2 

-1 -1 -1 Z4(8)- -1 26(8) 



- zl{8) - z|(8) - ^1(8) + 2zi (8)23(8)^4(8) 
1 

v/1 - zl{8) - z|(8) ~ z|(8) + 222(8)25(8)^6(8) 
1 1 1 

T~b^KJ T^~B^ T^M^ ' 



(23) 



/9(ijfc/) = (24) 



ho{ijklm) 



B,(.12)B^.(,34)s^.(,56) 

(?Tii+m2)^ 

('•"3 + '"4)^Aj ('"5 + '"6)^Afc 



rf44 g'(^o')P.(44) / ^46 Gli4')Pk{s'5e) 

34 — ■'0 J *56 ~ ^0^ 

+ 1 +1 +1 +1 Z5(10) + 

1 f dzi(lO) f dz2(10) /■ dz3(10) /■ dz4(10) 



X 



-1 -1 -1 -1 25(1-)- 

1 

- Z^(IO) - ziilO) - Zim + 221(10)Z4(10)Z5(10) 

1 1 



l-i?((s'23)l-5mK5) 



(25) 



where i, j, k, I, m correspond to the diquarks with the spin-parity = 0"*", 1"^ and mesons with the spin-parity 

jp = o-,i-. 

The other choices of point sq do not change essentiahy the contributions of a^, therefore we omit the indices s^j'. 
Since the vertex functions depend only shghtly on energy, it is possible to treat them as constants in our approximation. 

We can pass from the integration over cosines of the angles to the integration over the subenergies [2l[. In the 
relativistic invariant solution the center of mass of particles 1, 2 by the standard method is treated |21| . 

The system of graphical equations Fig. 1 is determined by the subamplitudes Ai, A2, A3. BB states {A2) are 
constructed without the mixing of quarks and antiquarks. Therefore we did not use the three mesons and meson plus 
tetraquark states. But the subamplitudes Ai and A3 contain the quark-antiquark pairs. Then the algebraic equations 

(13) - (17) take into account the contributions reduced amplitudes a\ , a\ , a\ , a\ ^ , a\ ^ ^ . The solution 
of the system of equations are considered as: 



a,{s)^F,{s,\,)/D{s), (26) 

where zeros of D{s) determinants define the masses of bound states of baryonia. 

We have analyzed in the subamplitudes of a quark and an antiquark and did not obtain the bound state with the 
model parameters. 

As example, we consider the equations for the quark content uuuddd with the isospin / = 3 and the spin-parity 
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= 3~ (Fig. 1). The similar equations have been calculated for the isospin / = 0, 1, 2, 3 and the spin-parity 
= 0^, 1^, 2^, 3^. We take into account the u and d quarks. 
The functions /i, I2, I3, In, Iq, Is, Iq, Iio determine the interaction of the quarks and the antiquarks. These 
functions take into account the contributions of the Chew-Mandelstam functions, which are constructed in the model 
for the quark-antiquark pairs with various quantum numders using the unitarity condition. 



III. CALCULATION RESULTS. 



The poles of the reduced amplitudes ai correspond to the bound states and determine the masses of the baryonia. 
The dynamical mixing between the subamplitudes of baryonia is considered. We derived the relativistic six-quark 
equations in the framework of the dispersion technique. The pair quarks amplitudes qq — > q g and qq — )■ qq are 
calculated with the dispersion N/D method using the input four-fermion interaction |22l . |23| with the quantum 
numbers of the gluon [24, 25]. 

The model under consideration proceeds from the assumption that the confinement radius is sufficiently larger than 
constituent quark radii as well as the radii of the forces which bound the low-lying hadrons. It means that quark 
interaction forces are the two component ones. The long-range component is due to the confinement. In the present 
paper, when the low-lying hadrons are considered, the long-range component of the forces is neglected. 

We manage with the quarks as with real particles. However, in the soft region, the quark diagrams should be treated 
as spectral integrals over quark masses with the spectral density p(m^): the integration over quark masses in the 
amplitudes puts away the quark singularities and introduces the hadron ones. We can believe that the approximation 
pijn?) 5[m? — m^) could be possible for the low-lying hadrons. We hope that this approach is sufficiently good for 
the calculation of the low-lying baryonia being carried out here. The problem of the distribution over quark masses 
is important when one considers the high-excited states. 

The four-quark interaction is considered as an input: 



gv mn^^qf ■ (27) 



Here // is the unity matrix in the fiavor space (u, d). A are the color Gell-Mann matrices. 
We introduce the scale of the dimensional parameters [25| : 



4A(ifc) 
TT- (mi + mky 



g=—gy, A=—^-±—. (28) 



Here and mj; are the quark masses in the intermediate state of the quark loop. Dimensionless parameters g and 
A are supposed to be the constants which are independent of the quark interaction type. In the case under question 
the interacting pairs of particles do not form the bound states. The attraction of quark-antiquark and quark-quark 
pairs is not enough for the construction of the bound state. This is similar to the case of the four-quark systems 32]. 
Therefore, the integration in the dispersion integral run from Am? to A. 

The quark masses of model rriu^d = ilOMeV coincide with ordinary baryon ones '21I]. The model in question has 
only two parameters: the cutoff parameter A = 11 and the gluon coupling constant g — 0.314. These parameters are 
similar to the previous paper [21| ones. 

The estimation of theoretical error on the baryonia masses is 1 MeV. This result was obtained by the choice of 
model parameters. 

Some authors had investigated the tetraquark system with the three-body qqq and qqq interaction, whose existence 
has no direct effect on the ordinary hadron states [33.., 34]. But assuming the most general two-body quark Hamiltonian 
[35I l36j . Pirjol and Schat derive universal correlations among masses and mixing angles which are valid in any model 
for quark interactions containing only two-body interactions. Deviation from these predictions provide no evidence 
for the presence of spin-flavor dependent three-body quark interactions. 

Our model is based on the three principles of unitarity, analyticity and crossing symmetry. The principle of unitarity 
are applied to the two-body subenergy channels. 

Further experimental and theoretical efforts are required in order to satisfactory explain the presence of three-body 
quark interactions. 

In the Table I the calculated masses of nonstrange baryonia are shown. The contributions of subamplitudes to the 
six-quark amplitude are represented in the Tables II, HI, IV. The states (AA + Ah + nA + nh), (AA -I- Ap + pA+pp) 
and ( AA + An -I- pA + ph) with the isospins 7 = 0,1 and the spin-parity ~ 0^ possess the mass M ~ 1835 M eV . 
We predict the degeneracy of the some states. For the (AA +pA), (AA -j- An), (AA -|- nA) and (AA -f Ap) with 
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the isospins / = 1, 2 and the spin-parity = 0~ the mass M = 1928 MeV is obtained. The low- lying state uuuddd 
(AA) with the isospin / = 3 and the spin-parity = 1~ possesses the mass M — 1783 MeV . The dynamical mixing 
between the five subamplitudes (similar to the Fig. 1) is considered. Therefore the multiquark state will be stable. 
The states uuuuuu (AA) and dddddd (AA) with the isospin 7 = and the spin-parity = 0~ have the mass 
M = l9TiMeV. 

We predict the degeneracy of baryonia M{uudddd, I = 2) ~ M{uuuudd, I ~ 2) ~ M{uddddd, I = I) = 
M{uuuuud, I ~ 1). For the states M{uududd, I = 1) = M{uddudd, 7 = 0) = M{uuduud, I — 0) and 
M{uuuuuu, I — 0) — M{dddddd, 7 = 0) the degeneracy is also obtained. 

A somewhat simple picture of baryonium is that of a deuteron-like NN bound state or resonance, benefiting from 



the attractive potential mediated by the exchange of gluon [21[. We consider the influence of the contributions of 
quark-antiquark pairs. 

Entem and Fernandez, describing scattering data and mass shifts of pp levels in a constituent quark model, assign 
the threshold enhancement to final-state interaction [s^, HI] . Zou and Chiang find that final state interaction makes 
an important contribution to the pp near threshold enhancement [s^ . 

BES collaboration measured the mass X(1835) to be M = 1833.7Mey and its width to be F = 67.7 MeV [2]. 

The state with M = 1835 AleV is considered as pp state [4!] or the second radial excitation of rj' meson (40| . 

In our case this state have following content A A -I- Ap + pA + pp with isospin 7 = and spin-parity = 0~. 

We calculated the masses of baryonia with the isospin 7 = 0, 1, 2, 3 and spin-parity = 0~, 1~, 2~, 3~ (Table 

I). 
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TABLE I: S-wave baryonia masses. Parameters of model: cutoff A = 11.0, gluon coupling constant g = 0.314. Quark masses 
mu,d = AlQMeV. 



I 


Quark content (baryonia) 


J 


Mass (MeV) 
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1 
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rlrlrliliiii ( /\A\ 
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2290 


0; 1 


uddudd (A A + An + nA + nn) , 
uuduud (A A + Ap + pA + pp) ; 
uududd (AA + An + pA + pn) 


u 
i 

9 

3 


loOO 
1 7© /I 

1 S^t^l 
loOl 

2455 




uddddd (AA + nA), 





1928 


1; 2 


uuuuud (AA + Ap); 


1 
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1 770 
1 / / U 


uudddd (AA + pA), 


2 


1857 




uuuudd (AA + An) 


3 


2395 









2067 


3 


uuuddd (AA) 


1 
2 
3 


1783 
1938 
2290 



TABLE II: I J = 00, uuduud (1835 MeV"), A = 11.0, p = 0.314. 



Subamplitudes 


Contributions, percent 




3.7 


Ar 


3.7 


Ar 


9.0 


Ar 


7.8 


A^ 


7.8 


Af 


6.7 


AT' 


3.6 


Af' 


3.6 


Ar 


7.0 


Af 


6.8 


Ar 


6.8 


Af 


6.6 




2.4 


Ar'"" 


2.0 


Ar'"' 


2.6 




2.8 




4.0 




4.5 


3 


4.5 


^3 


4.3 




73.0 




9.8 


E^3 


17.2 



TABLE III: IJ = 33, uuuddd (2290 MeF), A = 11.0, g = 0.314. 



Subamplitudes 


Contributions, percent 


1 __ 


Q Q 


^1 


9.9 




25.4 




14.5 


^3 


40.3 




45.2 




14.5 


E^3 


40.3 



TABLE IV: IJ = 31, uuuddd (1783 MeV), A = 11.0, g = 0.314. 



Subamplitudes 


Contributions, percent 


Ar 


14.5 


Al"' 


14.5 


Ar 


34.8 


Ar 


25.2 




ii.l 


E^i 


88.9 




11.1 



TABLE V: Vertex functions and Chew-Mandelstam coefficients. 



i 


G1{sm) 




Pi 


.5i 


0^ diquark 
1^ diquark 
0^ meson 
1^ meson 


3 (3sfe() 
2s 


1 
2 

1 


1 (mfe-mi)^ 

2 (mfc+m;)^ 

imirjni 1 




1 (wife-m;)^ 


3 

3 (33fci) 
4g 


3 

1 
2 

1 


3(mfe+m;)^ 6 
1 (mh-mif 

imk-mi 1 


6 (rrifc+mj)^ 


1 ("ife-mj)^ 


3 


3 


3(mfe+m;)^ 6 


6 ("T-fc+'"i)^ 



